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We investigate a generalized two-dimensional Weyl Hamiltonian, which may describe the low- 
energy properties of mechanically deformed graphene and of the organic compound a— (BEDT- 
TTF)2l3 under pressure. The associated dispersion has generically the form of tilted anisotropic 
Dirac cones. The tilt arises due to next-nearest-neighbor hopping when the Dirac points, where the 
valence band touches the conduction band, do not coincide with crystallographic high-symmetry 
points within the first Brillouin zone. Within a semiclassical treatment, we describe the formation 
of Landau levels in a strong magnetic field, the relativistic form of which is reminiscent to that of 
graphene, with a renormalized Fermi velocity due to the tilt of the Dirac cones. These relativistic 
Landau levels, experimentally accessible via spectroscopy or even a quantum Hall effect measure- 
ment, may be used as a direct experimental verification of Dirac cones in a— (BEDT-TTF)2l3. 
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I. INTRODUCTION 

The discovery of a particular quantum Hall effect in 
grapheneii^ has shown that the low-energy electronic 
properties in this two-dimensional (2D) carbon crystal 
are described not in terms of a Schrodinger-type wave 
equation but by a relativistic Dirac equation^ Due to 
a 7t-band, which shrinks at half-filling to two inequiva- 
lent points at the corners of the first Brillouin zone (BZ), 
the electronic energy dispersion is almost linear result- 
ing in Dirac cones. This is reminiscent of the case of 
massless relativistic particles, where the speed of light c 
is replaced by a Fermi velocity vf, which is roughly 300 
times smaller than c. 

Another material where Dirac cones are expected 
to occur is the organic 2D compound a— (BEDT- 
TTF)2l-3 under pressure4^£ The relativistic behavior 
of the carriers may be at the origin^ of an experimen- 
tally observed T 2 dependence of the carrier density. 7,8 
Whereas in graphene, the Dirac cones at the corners 
of the first BZ are isotropic, they are situated within 
the first BZ in a-(BEDT-TTF) 2 I 3 , strongly anisotropic, 
and tilted in the wave- vector energy space (k, E) The 
electronic properties are described by a generalized Weyl 
Hamiltonian with terms linear in the 2D wave vector k. 
However, in contrast to graphene, there is yet no direct 
experimental evidence for the presence of Dirac cones in 
a— (BEDT-TTF)2l3 or whether the system is simply a 
narrow-gap semiconductor. 

In the present paper, we study the structure of the 
generalized Weyl Hamiltonian, which yields energy dis- 
persions in form of tilted anisotropic Dirac cones. In 
the presence of a strong magnetic field, the dispersion 
is quantized in relativistic Landau levels (LLs), with 
the characteristic ±\/nB behavior known from graphene. 
The tilt and the anisotropy of the Dirac cones give rise 
to a renormalization of the effective Fermi velocity and 
therefore of the typical LL spacing. 

One example of a 2D system described by such gen- 
eralized Weyl equation may be the above-mentioned or- 



ganic material a— (BEDT-TTF) 2 l3. We show, within an 
effective tight-binding model on an anisotropic triangu- 
lar lattice with two atoms per unit cellj2 that the tilt- 
ing of the Dirac cones is due to next-nearest-neighbor 
(nnn) hopping, which may be in a— (BEDT-TTF) 2 l3 on 
the same order of magnitude as nearest-neighbor (nn) 
hopping i 10 ' 11 A necessary condition for nnn hopping to 
cause a tilt of the Dirac cones is that they are situated at 
points in the first BZ different from those of high crys- 
tallographic symmetry, such as its corners. Furthermore, 
we show that it may equally apply to graphene when 
the Dirac points, D and D' move away from the high- 
symmetry points K and K' at the corners of the first 
BZ. In this case the wave- vector expansion of the nnn 
term yields a linear contribution, whereas it is quadratic 
when the Dirac points coincide with the BZ corners K 
and K'. Such motion of the Dirac points may indeed 
be induced by a quinoid-type lattice distortion^ 2 - of the 
graphene sheet. However, we show that the tilt of the 
Dirac cones is much less pronounced than in a— (BEDT- 
TTF) 2 l3. Alternatively, this motion of Dirac points may 
be studied in cold atoms in an optical lattice where one 
may deform the honeycomb lattice and fine-tune the nn 
and nnn hopping parameters with the help of the laser 
intensities, wavelengths, and relative orientation. 13 

The paper is organized as follows. We start with a the- 
oretical discussion of the generalized Weyl Hamiltonian 
in Sec. II. Sec. Ill is devoted to the LL formation in a 
strong magnetic field, for the case of tilted Dirac cones. 
Possible experimental realizations in distorted graphene 
and a— (BEDT-TTF) 2 I 3 are discussed in Sec. IV, which 
we conclude with an analysis of a possible quantum Hall 
effect in a-(BEDT-TTF) 2 I 3 . 



II. GENERALIZED WEYL HAMILTONIAN 

We consider a model of two-spinor fermions restricted 
to a 2D space. Whereas the two-spinor form is in gen- 
eral dictated by relativistic invariance in two space di- 
mensions, it naturally arises in the condensed matter 




situation of a lattice with two inequivalent sites. The 
most general Hamiltonian linear in the 2D wave vector 
k = (k x ,k y ), is given by the "generalized Weyl Hamilto- 



H = 



M=0,...,3 



• ko^ 



(1) 



in terms of the velocities v 



„ - and the 2x 2 Pauli 

matrices u u = l,a = (a 1 , a 2 , a 3 ). Here and in the fol- 
lowing parts, we choose a unit system with h = 1. Both 



2D space components of the velocities, 



(u£,uf,vf,'if) and v» = {v%,v") = [vl,v\,v%,vl), are 
in themselves vectors in the 4D spin space [the space 
of SU(2) matrices] spanned by the Pauli matrices. The 
usual 2D Weyl Hamiltonian, which describes for instance 
low-energy massless electrons in graphene^ is included 
in (Q} if one considers v = v 4 = 0, vi = (vp,0), and 
V2 = (0, vp), in terms of the Fermi velocity Vf- 

Although, at first sight, the Weyl Hamiltonian is de- 
scribed by eight different parameters, given by the four 
two-component velocities v^, it is indeed overspecified. 
In order to illustrate this point, we rewrite the Hamilto- 
nian (TT]) in a different manner, 



H = v -kcr° + {v x k x 



(2) 



One may get rid of two parameters (V3 = 0) by choosing 
the 3-quantization axis in the SU(2) space perpendicular 
to the vectors v x and v 1 * ' . 

This point is indeed remarkable and needs to be dis- 
cussed in the light of graphene physics. In this case, a 
constant tr 3 term breaks the inversion symmetry of the 
honeycomb lattice, e.g. due to a different on-site energy 
of the two triangular sublattices. Usually, this gives rise 
to a mass term and breaks the particle-hole symmetry. 
In the generalized Weyl Hamiltonian, this is not the case 
because the a 3 term is linear in the wave vector and 
therefore does not affect the zero-energy state at k = 0. 

One may furthermore reduce the number of relevant 
model parameters by a simple rotation of the 2D frame 
of reference, accompanied by a unitary transformation 
in the SU(2) space, which leaves the 3-quantization axis 
invariant. One, thus, obtains the "minimal" Weyl Hamil- 
tonian 



H = w • qcr° 



w x q x a 



v v q v o» 



(3) 



in terms of the four effective velocities Wo = {wo x ,wo y ), 
w x and w y . A detailed discussion of the involved trans- 
formations and a derivation of the exact expressions for 
the effective velocities may be found in the Appendix [A"l 
The diagonalization of the minimal Weyl Hamiltonian 
yields the energy dispersions 



e A (q) = w ■ q + A 



wWx 



2 y q 2 y , 



(4) 



where A = ± plays the role of the band index. 



For w x 



w„ 



Vp and w 



isotropic model, which applies e.g. 



0, one obtains the 
to the low-energy 




FIG. 1: Energy dispersion Q for the special choice of w x = 
w y — 1, and wo = (0,0.6), in natural units. The Dirac cone 
is tilted in the y-direction. 



electronic properties in graphene: the Fermi velocities 
are the same in the x- and y-direction. The rotational 
symmetry is broken if w x ^ w v (anisotropic model) . Such 
case may be obtained e.g. if the graphene sheet is con- 
strained by a uniaxial pressure, as is discussed in Sec. IV 
A. For wo 0, the Dirac cones are tilted away from the 
z-axis, as is shown in Fig. [TJ 

Notice that not all values of the tilt parameter wo are 
indeed physical. In order to be able to associate A = + 
to a positive and A = — to a negative energy state, one 
obtains the condition 



(^Oy) <L 



(5) 



Unless this condition is satisfied, the iso-energetic lines 
are no longer ellipses but hyperbolas. Notice that, here, 
we aim to use the generalized Weyl Hamiltonian ([I]) and 
its resulting energy dispersion ^ to describe the elec- 
tronic properties of particular 2D materials. Although it 
may be interesting to speculate about the resulting prop- 
erties of a model that violates the condition ([5]), we are 
not aware of any physical example which might corre- 
spond to such a case. 

In a 2D lattice system with valley degeneracy, a gen- 
eralized Weyl Hamiltonian may describe the low-energy 
excitations in different valleys separately. In the remain- 
der of this paper, we will in general only consider a single 
valley (and explicitly mention the inclusion of the twofold 
valley degeneracy when needed). Note also that we do 
not consider the true electron spin and do not include 
the corresponding twofold spin degeneracy. 

In order to discuss the symmetry properties of the gen- 
eralized Weyl Hamiltonian ([3]), it is convenient to intro- 
duce the unitary and Hermitian chirality operator 



C = 



w 2 x q 2 x 



w 2 y q 2 y 



(6) 



3 



which commutes naturally with the Hamiltonian. The 
associated eigenvalues are a = ±1 and coincide with the 
band indices a — A. As exemplified in Sec. IV, this is 
generally not the case in a physical condensed-matter sit- 
uation - the Weyl Hamiltonian corresponds to the effec- 
tive model at Dirac points, where the conduction band 
touches the valence band; these Dirac points occur in 
pairs, at inequivalent points in the first BZ, which yields 
a twofold valley degeneracy. In this case, the effective 
model is rather given by £H, where £ = ± denotes the 
two valleys, and the relation between band index, chiral- 
ity, and valley index is given by 



(7) 



In the present discussion, we may however identify the 
band index with the chirality, for simplicity. 
The eigenstates of the chirality operator are 



1 / e -i<Pk 



V2\ 



a 



(8) 



where tan 0k = w y k y /w x k x . These eigenstates are also 
the natural eigenstates for the generalized Weyl Hamil- 
tonian. 



III. TILTED DIRAC CONES IN A MAGNETIC 
FIELD 

We use the Peierls substitution to obtain the general- 
ized Weyl Hamiltonian in a magnetic field 



q — ► II = q + eA, 



(9) 



where A is the vector potential that generates the (uni- 
form) magnetic field Be z = V x A perpendicular to the 
2D plane. With the help of the ladder operators 



at = 



y/2w x Wy 

h 

^2w x Wy 



{W X U X - IWyXly) , 

{w x Il x + iWyTLy) , (10) 



in terms of the magnetic length Is = 1/ V eB, one obtains 
the Hamiltonian 



H B = 



^2w x w y ( ^(ae i( c + aU'^) 



where we have defined 

j» - w ** | . : w Qy 



2f (oe* +at e -^) 
(11) 



w e 



in terms of the effective tilt parameter 



w 




WOyY 

Wy J 



(12) 



Instead of the full solution of the Hamiltonian (fTTj) , we 
consider the effect of the magnetic field in a semiclassical 
treatment. The Onsager relation 14 states that the sur- 
face S(e) enclosed by a trajectory of constant energy e in 
reciprocal space is quantized as 



(2tt) 2 / de'p(e') = 27r(n + 7 ), 



S(e)l% 



where n is an integer denoting the energy level which co- 
incides with the Landau level in the full quantum treat- 
ment. The additional contribution 7 is related to a Berry 
phase acquired by an electron during its cyclotron orbit. 
Usually, one has 7=1/2 except if there is an extra Berry 
phase of 7r, which in our case yields 7 = 0, as in the case 
of graphene with no tiltj 1 ^ If one considers a density of 
states which scales as p(e) cx e a , the energy levels thus 
scale as 



[B(n+7)] 1/(1+a) , 



(13) 



in the large- n limit. In usual (non-relativistic) 2D elec- 
tron systems, one finds a constant density of states, i.e. 
a = 0, and 7 = 1/2. The scaling of the conventional 
Landau levels is therefore e n cx B(n + 1/2). In the rela- 
tivistic case of electrons in graphene, the density of states 
vanishes linearly at the Dirac points, and one therefore 
obtains e„ cx V Bn because a — 1 and 7 = 0. The rela- 
tion (fT3)l has been generalized to the case of a spatially 
anisotropic density of states by Dietl et al*£ 

From the scaling argument (|13|) in the large-n limit, 
one may notice that the -B-field scaling of the levels must 
be the same as the n scaling. Furthermore, one sees from 
the quantum Hamiltonian (jlip that the energy must scale 
as 1/Ib cx \f~B. Therefore, the energy levels must obey, 
in the large-n limit, the equation 



(14) 



as in the case of the Weyl equation for massless charged 
particles, such as in graphene, apart from a renormaliza- 
tion of the Fermi velocity. 

The renormalization of the Fermi velocity may be ob- 
tained from the calculation of the density of states. The 
total number of states below a given energy e within the 
positive energy cone is given by 



N+(e) 



1 



(2Tr) 2 w x w y J e+i g) 



dq x dq y 



2nv* F 2 2 ' 



where we have defined q x / y = w x / y q x / y , and the renor- 
malized Fermi velocity is written in integral form, 



1 

TP? 



1 



2tt 



1 



w x w v Jq 2n (1 + wo cos ( 



(15) 



in terms of the effective tilt parameter (TT21 . One notices 
from Eq. (TT5|) that if the condition (0, |wq| 2 < 1, is 
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FIG. 2: Renormalized Fermi velocity v F / ^/w x w y 
as a function of the effective tilt parameter w = 
■J (wox/w x ) 2 + (wQy/vjy) 2 . The Fermi velocity vanishes 
for w — 1, where the orbits change from ellipses to 
hyperbolas. 

not satisfied, the expression under the integral diverges 
because the denominator may become zero. This result 
is not suprising because the Onsager quantization rela- 
tion, which yields the energy levels (fT4")l is only valid for 
closed orbits, given e.g. by the elliptic isoenergetic lines. 
As already mentioned, the orbits for \wq\ > 1 are open 
hyperbolas, and the expression (|14[) is no longer valid. 

The density of states is obtained by differentiation of 
the number of states, 



which is the concise expression for both the positive and 
negative parts of the tilted Dirac cones. 

The \JnB behavior of Eq. ([XI)) is, strictly speaking, 
valid only in the large- n limit. However, usually it yields 
extremely good estimates for the levels down to values as 
small asn= 1. Special care is needed for the discussion 
of the n = level, which requires a quantum treatment 
of the Hamiltonian (fTT|) . In the following, we discuss the 
fate of the zero-energy Landau level. 

The behavior of this level may be understood with the 
help of the quantum treatment of the Hamiltonian for 
wo = 0. In this case, the expression (|14|) is exact with 
Vp = ^/w x Wy, which is also the Wo = 0-limit of the 
expression (fTSj) . There exists thus a zero-energy level 
for n = 0, which has the same degeneracy, Nb as all 
other levels (A, n), in terms of the number of flux quanta 
Nb = AB/(h/e) threading the total surface A. 

For non-zero values of Wo, the Hamiltonian (|11[) may 
not be diagonalized by a simple canonical transforma- 
tion. However, the Hamiltonian (fTTj) is transformed as 
Hb — * — Hb under space inversion, r — * — r, as shown 
in the Appendix [Bj This implies that the energy spec- 
trum is symmetric around zero energy (see Appendix IB"]) . 
Therefore, starting from wo = and adiabatically turn- 
ing on wo 7^ 0, there are only two possibilities for the 
evolution of the zero-energy level: (i) either it remains at 



zero energy or (ii) it splits into (at least) two sublevels + 
and _ , which are symmetric around zero energy. How- 
ever, splitting of the zero-energy level into sublevels can 
be excluded on account of the degeneracy of this level. 
Indeed, when Wo = 0, the exact degeneracy of the zero 
energy (n = 0) Landau level is given by Nb (remem- 
ber that we only consider a single valley here). When 
wo 7^ 0, it can, therefore, not split since this would indeed 
lead to an unphysical doubling of the number of quantum 
states because each level, + and 0~, would have to be 
Nb times degenerate. Therefore, for all magnetic field 
strength, there exist a zero-energy Landau level. The ex- 
plicit expressions for the zero-energy wave functions may 
be found in the Appendix [Cl Notice, that this is consis- 
tent with the semiclassical spectrum with 7 = 0. 

In the above treatment, we only considered a single 
valley. We note however that the magnetic field might 
introduce a coupling between the two valleys. In such 
a case, we do not exclude a parity anomaly which con- 
sists of a different behavior of the n = level at the two 
inequivalent Dirac points at non-zero wave vectors in a 
lattice model. In this case, space inversion would involve 
the low-energy Hamiltonians at both Dirac points, and 
the spectrum is only symmetric around zero energy if one 
accounts for both valleys. The parity anomaly is, how- 
ever, expected to play no physical role in the continuum 
limit with a/Is —* 0, where a is the lattice spacing. 

In conclusion, we have obtained the semiclassical spec- 
trum of Landau levels [see Eqs (fT4"|) - (fT2] , valid when 
n ^> 1) and checked that the zero energy level (n = 0) 
indeed exists in a full quantum treatment. Based on this 
two calculations, we expect the semiclassical spectrum to 
be a very good approximation to the true quantum spec- 
trum of Landau levels, for all n. This is one of the main 
result of the present paper. 



IV. PHYSICAL EXAMPLES OF TILTED DIRAC 
CONES 

After this rather technical discussion of the general- 
ized Weyl Hamiltonian and tilted Dirac cones, we dis- 
cuss, here, two physical systems which may display these 
properties. We find that whereas the tilt of the Dirac 
cones is well pronounced and thus strongly affects the 
Landau level quantization in a— (BEDT-TTF^Lj, it is 
much more difficult to induce a tilt in graphene via a lat- 
tice deformation. However, a quinoid-type lattice defor- 
mation is also discussed for pedagogical reasons because 
the general physical origin of tilted Dirac cones becomes 
transparent. 



A. Quinoid-type graphene under uniaxial strain 

As a first example, we consider a graphene sheet which 
is deformed in one of its principle symmetry axes. This 
particular deformation results in a quinoid variety of the 
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FIG. 3: Quinoid-type deformation of the honeycomb lattice 
- the bonds parallel to the deformation axis (double arrow) 
are modified. The shaded region indicates the unit cell of 
the oblique lattice, spanned by the lattice vectors ai and a2. 
Dashed and dashed-dotted lines indicate next-nearest neig- 
bors, with characteristic hopping integrals t nnn and t' nnn , re- 
spectively, which are different due to the lattice deformation. 



honeycomb latticed We treat its electronic properties 
within the tight-binding approximation. Starting from 
the graphene honeycomb lattice, with equal bond length 
a ~ 0.14 nm and equal nn hopping energy t ~ 3 eV, 
the bond length and hopping energy are modified in the 
deformation axis (see Fig. [3J, 



8a 



and 



t 



8t r 
t' = t + —Sa, 
da 



and kept unchanged otherwise. We call e = 8a/ a the 
relative strain. Here, we consider a moderate deforma- 
tion, |e| <C 1, such that one may linearize the hopping 
energy around its nondeformed value t, and dt/da ~ —5 
eV/A^LiS. This value agrees with an evaluation based on 
Harrison's law-i^ according to which t = Ch 2 /ma 2 , where 
C is a numerical prefactor of order one. Derivation with 
respect to a yields 



-4.3 eV/ A. 



(17) 



dt 2t 
da a 

For simplicity and as a first approximation, one may keep 
the bond angles fixed at 27r/3. The underlying Bravais 
lattice is no longer triangular but oblique with the basis 
vectors 



ai 



\f?>ae x 



and 



a 2 



V3 



3 

—a + 8a 



and the reciprocal lattice is spanned by the vectors 

47re„ 



2tt 



y/Za 3a + 28a 



and a5 



3a + 28a 



Furthermore, we take into account nnn hopping, with a 
characteristic energy oS^ t nnn ~ O.li in the undeformed 
horizontal axes. The deformation yields, in the same 
manner as for the nn hopping energies, different hopping 
energies for the other directions (see Fig. [3]), 



in 



dt n 



Oa 



l -8a. 



The tight-binding model may be described by the 
Hamiltonian 



ff=£(«U)w q (£;) 



(18) 



in reciprocal space, where a q ^ and are the Fourier 
components of the annihilation (creation) operators on 
the A and B sublattices, respectively. The Hamiltonian 
2x2 matrix 



W q = 



h'(q) h*(q) 
h(q) h'(q) 



is given in terms of the elements 



h(q) 



-t 



D i(q y + V3q x )a/2 , e i(q y -s/3q x )a/2 



—t'e~ iqy ( a +^ a ) 

q y a V3q x a 
= — 2icos-|-cos — t cos [q y {a + 8a)\ 



(19) 



■ Jo,- 1y a V%q x a , . . . 
-i < 2rsm -|- cos — t sin [q y (a + 8a)\ 



and 



h'(q) 



2t nnn cos V3q x a 

V3q x a 



-^t'nnn { COS 



2 + % ( 2° + 6a 



COS 



V3q x a /3 
2 — + 1y [ 2 a + Sa 



(20) 



The energy dispersion is obtained from the eigenvalues 
of H c 



e A (q)=/i / (q)+A|/ l (q)| 



(21) 



and is plotted in Fig. [4] for a deformation of 8a/ a = 0.4. 
The two bands, A = + and A = — , touch each other at the 
Dirac points q D , which are obtained from the condition 
h(q D ) = 0^ 







and 



9x a 



£ — = arccos 
V3 



(22) 



where £ = ± denotes the two inequivalent Dirac points 
D and D' , respectively. In the absence of any distortion, 
the Dirac points D and D' coincide with the crystallo- 
graphic points K and K' , respectively, at the corners of 
the first BZ. The distortion makes both pairs of points 
move in the same direction due to the negative value of 
dt/da. However, unless the parameters are fine-tuned, 
this motion is different, and the two pairs of points no 
longer coincided 

The low-energy properties of electrons in a quinoid- 
type distorted graphene sheet are described by the lin- 
earized model around the Dirac points, which is exactly 
of the form ([3]) of the Weyl Hamiltonian, 



= C ( w o • k cr° + w x k x a x + w y k y a v ) 



(23) 
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FIG. 4: Energy dispersion of the quinoid-type deformed the 
honeycomb lattice, for a lattice distortion of 8a/a = 0.4, with 
t = 3 eV, trmn/t = 0.1, dt/da — —5 eV/A, and dt nnn /da — 
—0.7 eV/A The inset shows a zoom on one of the Dirac points, 
D' . 



parabolic correction, which breaks the original electron- 
hole symmetry£i22 

To summarize, in order to obtain tilted Dirac cones in 
graphene, two ingredients are required: (i) nnn hopping, 
which generates the diagonal components h'(q) in the 
Hamiltonian (fT5)) ; and (ii) for a linear contribution aris- 
ing from this term, the Dirac points D and D' need to be 
shifted away from the high-symmetry points K and K' . 
This shift may be obtained by constraining the graphene 
sheet into such a quinoid type. 

In the presence of a magnetic field, the LL spacing is 
affected by the deformation because the Fermi velocity is 
renormalized according to Eq. (|15l) . 



1 



(25) 



for small values of the effective tilt parameter wq. It may 
be evaluated from the model parameters, 



w 



tnnn. Sm 29 tL 



t 



t 2 



("lira 



sin 9 

- t% 



f 



(26) 



with the effective velocities 



= V3ta sin t 



Wy 

w 0x 

W y 



25a 



-t'a ( 1 , 
2 V 3a 



2-\/3 (t nnn a sin 29 + t' nnn a sin ( 
0, 



(24) 



where we have defined 9 = arccos(— t'/2t). The corre- 
sponding energy dispersion is independent of £ , which is 
at the origin of the twofold valley degeneracy. In order to 
obtain the concise form of Eq. (f2"5| . we have chosen the 
spinor representation (jPa^b) at the £ = + Dirac point 
and (iPb^a) for £ = — , i.e. interchanged the sublattice 
components at D' . As mentioned in Sed. II, the relation 
between the band index A, chirality a, and valley index 
£ is given by Eq. (JT]), A = £a, due to the global sign £ in 
the Hamiltonian (j2li|) . The constant term ft/(q = ^q )! 
has been absorbed in a renormalization of the chemical 
potential, the position of which is determined by the elec- 
tronic half- filling of the graphene sheet. 

One notices from the Eqs. (|2^|) that the quinoid-type 
distortion yields an anisotropy in the Fermi velocities, 
w x 7^ w y , and that the Dirac cones are tilted due to 
wq x 7^ 0. The isotropic graphene model is retrieved at 
5a = - one has then w x — w y = vf = ita/2 ~ 6.3 
eVA and w 0x = w 0y = because t = t', t nnn = t' nnn , 
and sin# = \/3/2 = —sin 29, in the undeformed case. 
Without deformation, nnn hopping therefore does not 
affect the energy dispersion at linear order, but only at 
second order. This is due to the fact that the Dirac 
points are then situated at the high-symmetry crys- 
tallographic points K and K' . Indeed, this yields a 



In order to estimate t' nnn , we use the "atomic orbitals 
overlap law" familiar in the context of the extended 
Hiickel model, 23 



t nnn {b, a) ps t{a)e 



-(b-a)/d(a) 



where a is the nn distance, b is the nnn distance, and 
d w a/3.5 w 0.4 A is a caracteristic distance related 
to the overlap of atomic orbitals. In the undeformed 
graphene b = a\/3, whereas in the quinoid type graphene 
b' = 6(1 + e/2) and a' = a(l + e). This gives t' nnn = 
t nnn (l -2e + be/2d) and = t(l - 2s). Therefore, the 
effective tilt parameter is given by 



b t nn 



Wo S3 -- 

a 



t 



0.6e. 



As the correction to the Fermi velocity appears as 1 — 
3?ig/4 [see Eq. ([2"5])]. this effect remains extremely small, 
and the tilt affects the LL spacing in a negligible manner. 

The main contribution to the renormalized Fermi ve- 
locity therefore arises not from the tilt of the Dirac cones 
(effect of order e 2 ), but from the anisotropy in the Fermi 
velocities (effect of order e) , and one finds 



VF 



1 



dt 5a 
da t 



vf 



(l-§), (27) 



which may yield an experimentally observable effect in 
the percent range for a strain of e ~ 10%. 

/From an experimental point of view, such quinoid- 
type deformation may be realized if one uses a piezo- 
electric substrate, on which the graphene sheet is posed, 
instead of the most commonly used SiC>2. Another 
possibility would be to use a mechanical deformation 
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of the underlying substrate. Such bending has been 
exploited e.g. to investigate carbon nanotubes under 
strain.— More recently, graphene on polydimethylsilox- 
ane (PDMS) has been put under uniaxial strain by bend- 
ing of the PDMS. 2 » The elastic regime in graphene re- 
quires that the strain is smaller than 10% and the rupture 
occurs around 20%. Therefore an upper bound for e is 
certainly 10%. 



B. Organic 2D compounds 

Another example of a 2D metal, where tilted Dirac 
cones may occur, is the layered organic compound 
a-(BEDT-TTF) 2 I 3 under (uniaxial) pressurei^ Each 
layer may be described by an oblique lattice with four 
sites per unit cell, and the electronic filling is 3/4. In 
the vicinity of the Fermi energy, only two out of the four 
bands are relevant for the low-energy electronic prop- 
erties. It has indeed been shown that the band struc- 
ture may be modeled with great precision within a tight- 
binding model on a half-filled anisotropic triangular lat- 
tice with nn and nnn hopping, where each site corre- 
sponds to a dimer.— This is a natural assumption for k- 
and A-(BEDT-TTF) 2 l3, where there exists one hopping 
energy which is largely enhanced with respect to the oth- 
ers. In contrast to these compounds, the assumption may 
seem hasardous at first sight in the case of a— (BEDT- 
TTF)2l3, where there is no such clearly enhanced hop- 
ping energy, such that the dimerization is expected to be 
rather weak. Furthermore, these organic materials ex- 
hibit strong electronic correlations, and a tight-binding 
calculation for quasi-free electrons sweeps a lot of in- 
teresting physics under the carpet. However, the high- 
pressure limit corresponds to a regime where the elec- 
trons are less strongly correlated and where interaction 
effects may be taken into account via renormalized effec- 
tive hopping parameters £ 

The tight-binding model on the anisotropic triangular 
lattice is depicted in Fig. [5] The nn are situated at the 
vectors ±ai and iag, with 



Si = ^(e x +e y ) 



and <S 2 = ^{e x - e y ), 



which connect sites on the different sublattices, A and B 
and the vectors 



ai = Si + S? 



and 



a 2 



Si - S 2 



span the underlying Bravais lattice, which is chosen to 
be a square lattice, for simplicity. Notice that the lattice 
may also be viewed as an anisotropic 2D NaCl lattice 
(two inequivalent interpenetrating square lattices). The 
bond length is set to unity, a = 1. The nn hopping 
energies are ti and t\ in the directions ±aj., and t 2 and 
t' 2 in the directions z fa 2 , respectively. The nnn hopping 
energy is t nnn . 

The effective tight-binding model may be written in 
the same manner (|18p as for the case of quinoid-type 




FIG. 5: Anisotropic triangular lattice model, with four differ- 
ent nn hopping energies, ti, t'x, t<z, and t' 2 and the nnn hopping 
energy t nnn . The unit cell with two inequivalent sites is rep- 
resented by the shaded region. The sites of the A and B 
sublattices are depicted by the filled and open circles, respec- 
tively. 



graphene, with the matrix elements 



h(q) = 2 



{ti + t[) cos 



Qx + q v 



+2i 



(ti — t[) sin 



2 



(t 2 + 1' 2 ) cos 



q x - q y 



+ (*2 — t'2) sm ' 



and 



h'(q) = 2t nnn cosq y . 



The energy dispersion is obtained from Eq. (|21[) . and the 
position of the Dirac points is calculated from 



2 Qx 
tan 2 — 

2 



-,D 



tan 



(t[-t 2 f-(ti-t' 2 ) 2 
(t' 1+ t' 2 ) 2 -(t 1+ t 2 f 
{t' 1 -t' 2 f-{t 1 -t 2 f 



One may directly see that the r.h.s of both equations 
must be positive in order to have a pair of Dirac points 
(q and -q D ) within the first BZ, —n/2 < q x ,q y < n/2. 

An expansion around the Dirac points yields the gen- 
eralized Weyl Hamiltonian (fTJ), 

2 
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in terms of the velocities 



v Qx = 0, v 0y = -2t nnn sin q: t 



< = ( 


t'l+h) 


i sin 




+ 

2 


9y +( 


4+^) 


. 9? 
sin 


-1? 

2 


v\ = ( 


t'x + ti) 


i sin 


a? 


+ 
2 


<_ ( 




sin 


2 


V2 = ( 


t'i - h) 


i cos 


a? 


+ 
2 


<? + < 


4 - *a) 


cos 


-1? 
2 


vl = ( 


t'i - h) 


i cos 




+ 
2 




4 - fa) 


COS 


2 

(28) 



Here, we have used the same spinor representation as 
for quinoid-type graphene, i.e. we have interchanged the 
sublattice components when changing the valley. One 
notices that the Dirac cones are tilted only if the Dirac 
points are not situated at the border of the first BZ, 
qy = ir/2. This corresponds to the high-symmetry crys- 
tallographic points in graphene, and nnn hopping affects 
the effective model again only at second order in the ex- 
pansion around the Dirac points. 

The experimental evidence for (tilted) Dirac cones in 
a— (BEDT-TTF)2i3 compounds under pressure is yet 
rather weak. Whereas at ambiant pressure, the mate- 
rial is an insulator due to charge ordering, temperature- 
dependent transport measurements under high hydro- 
static pressure have revealed a T 2 dependence of the 
carrier density below 50 Kj££ as one would expect for 
relativistic electrons with a linear dispersion relationi^ 
It is, however, not clear whether the compound has, un- 
der these circumstances, a truely vanishing gap as for 
massless relativistic electrons or whether a tiny gap per- 
sists. Furthermore the T 2 dependence of the carrier den- 
sity is accompanied by a temperature-dependent mobil- 
ity, which results in an essentially constant conductivity 
over a large temperature ranged 

A more direct evidence for the relevance of Dirac 
cones in a— (BEDT-TTF)2i3 would be a measurement 
of the characteristic properties of relativistic quantum 
Hall physics, as in the case of graphene The follow- 
ing final part of this paper is devoted to the discussion of 
possible quantum Hall physics in a— (BEDT-TTF) 2 l3. 

Possible quantum Hall effect in a-(BEDT-TTF) 2 I 3 

Although it is a delicate issue to yield energy values 
for the hopping parameters t\, t' 11 t<i, t' 2 and t nnn from the 
overlap integrals in a-(BEDT-TTF) 2 l3r^ we expect 
that the good agreement between band-structure calcu- 
lations in the full model with four sites per unit cell and 
the anisotropic triangular lattice model^ yields the cor- 
rect orders of magnitude for the effective velocities (|28|l . 
Using the prescription proposed by Hotta^ and the over- 
lap integrals calculated by Mori et al.^ we may estimate 
ti = 36 meV, t' x = -86 meV, t 2 = -24 meV, t' 2 = -77 



meV, and t nnn = —60 meV. These values yield a pair of 
Dirac points at q° and ~q D , with q D = (2.52,-3.08), 
in units of the inverse lattice constant, which is on the 
order of 10 A.£Li±M with the help of Eqs. JSgJ, one 
thus obtains the effective velocities vf — —0.035 eVA, 
v\ = 0.315 eVA, uf = -0.222 eVA, v\ = -2.121 eVA, 
VQx = 0j and vq v = 0.074 eVA. One notices a variation by 
almost two orders of magnitude, and one may therefore 
expect rather large anisotropies. 

The effective velocities in the minimal model are cal- 
culated with the help of Eqs. (|A3|) . and one finds a 
rotation angle of 9 = 0.102 and the velocities w x — 
2.14 eVA, w y = 0.22 eVA, w 0x = -0.0075 eVA, and 
wo y = 0.736 eVA. The average Fermi velocity is there- 
fore y/w x w y = 0.69 eVA, which is roughly one order of 
magnitude smaller than that in graphene. The tilt pa- 
rameter (|12p is 

Wo = 0.33 

and thus much larger than in the case of a quinoid-type 
deformation of a graphene sheet. The tilt therefore leads 
to a reduction of the average Fermi velocity, and one finds 
from Eq. (]25p a renormalized velocity of 

v* F ~ 0.92^Tj~E^" ~ 0.63eVA. 

The renormalized Fermi velocity allows one to ex- 
tract the typical energy scale for the Landau levels in 
a-(BEDT-TTF) 2 I 3 , and one finds from Eq. (TJU) e A ,„ = 
XujQy/n, with a characteristic "cyclotron" frequency of 

uj* c = V2y^ ~ ZAy/B[T\ meV, (29) 

which is, due to the smaller Fermi velocity, roughly one 
order of magnitude smaller than that in graphene. How- 
ever, this energy scale is comparable to the cyclotron fre- 
quency in GaAs heterostructures {uc — 1.6S[T] meV), 
which are most commonly used in the study of quantum 
Hall physics. 27 One may therefore expect that a relativis- 
tic quantum Hall effectii 2 - could principally also occur in 
a-(BEDT-TTF) 2 I 3 if disorder does not prevent LL for- 
mation. 

Experimentally, thin (BEDT-TTF) 2 I 3 films have al- 
ready been synthesized. 28 Alternatively, one may hope 
that the exfoliation technique, 29 which has proven to be 
particularly successful in the fabrication of single-layer 
graphene sheets, also yields reasonably thin a— (BEDT- 
TTF) 2 I 3 samples. However, (BEDT-TTF) 2 I 3 crystals 
are generally of lower mechanical stability than carbon 
crystals, due to the relatively large lattice constants and 
the reduced binding energies. 

Apart from a direct measurement of a quantum Hall 
effect in a— (BEDT-TTF)2l3 compounds, one may probe 
the system via transmission spectroscopy in a magnetic 
field. This would allow for a direct measurement of the 
cyclotron frequency and for a check of the relativistic 
character of electrons in a-(BEDT-TTF) 2 I 3 . Trans- 
mission spectroscopy has indeed been successfully ap- 
plied to epitaxial^ and exfoliated 31 graphene and yields 
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a \fB{\Jn + 1 ± y/n) scaling of the transmission lines, 
as expected for the relativistic quantum Hall effect in 
graphene. 
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V. CONCLUSIONS 

In conclusion, we have investigated tilted Dirac cones 
in deformed graphene and the organic 2D material 
a— (BEDT-TTF)2l3- The low-energy electronic proper- 
ties are described by a generalized Weyl Hamiltonian, 
which may in both physical systems be derived from a 
tight-binding model on a lattice with two incquivalcnt 
sites. Whereas the presence of pairs of Dirac points is 
due to nn hopping, which couples neighboring sites on 
inequivalent sublattices, the tilt of the Dirac cones arises 
from nnn hopping if the Dirac points are shifted away 
from the points of high crystallographic symmetry in the 
first Brillouin zone. 

In the presence of a strong magnetic field, a semi- 
classical analysis yields the same structure of relativis- 
tic LLs as in non-deformed graphene, but with a renor- 
malized effective Fermi velocity due to the tilt and the 
anisotropy of the Dirac cones. Whereas this effect is ex- 
pected to be small in a quinoid-type deformation of the 
graphene, our estimates for the effective velocities for 
a-(BEDT-TTF) 2 I 3 indicate that the tilt yields a sig- 
nificant reduction of the effective Fermi velocity, which 
determines the LL spacing. The largest spacing of the 
— > +1 a nd —1^0 LL transitions is on the order of 
3Ay/B[T] meV, which is on the order of the (equidis- 
tant) LL spacing in GaAs heterostructures most com- 
monly used in quantum Hall effect measurements. Such 
measurements in a— (BEDT-TTF)2l3, as well as LL spec- 
troscopy, may be a possible experimental verification of 
the yet weakly corroborated presence of Dirac cones in 
a— (BEDT-TTF) 2 i3 . 



APPENDIX A: DERIVATION OF THE MINIMAL 
WEYL HAMILTONIAN 

In order to reduce the number of effective parameters 
in the Weyl Hamiltonian © , one rotates the 2D reference 
system in the physical space, 

k x = cos i9 q x + sin $ q y 

k y = — sin d q x + cos d q y , 
accompanied by a unitary transformation in the SU(2) 
space, 



U{6) 



cos - 1 
2 



i sin — a 
2 



which leaves the 3-quantization axis invariant and de- 
scribes a rotation in the rry-plane in the SU(2) spin space, 



cos a x + sin 6 a v 
— sin 9 a x + cos 9 a % 



(Al) 
(A2) 



If one chooses 



tan = 



and 



tan2i? 



vf sin $ + vf cos d 
v x sin •& + vK cos $ 



2(vfv\ + v$v%) 



1 

2v* 



>!) 2 
if 



K) 2 



K) 2 



m 2 ' 
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one obtains the "minimal" Weyl Hamiltonian ([3]). In 
terms of the original velocities, the minimal set of effec- 
tive parameters (the velocities Wo = (wq x , wq v ) , w x and 
w y ) reads 



w 0x 



Vq cos — Vq sin # , wo y = Vq sin $ + Vq cos i? 

K) 2 + K) 2 + K) 2 + K) 2 , 



K) 2 + K) 2 - K) 2 - K) 2 



+ (vfv\ + v% u|) 



(A3) 



K) 2 + {v* 2 f + (vf) 2 + K) 2 



{vf ) 2 + (vf) 2 - K) 2 - {vD 2 



+ (vfv\ + v x v\) 



APPENDIX B: TRANSFORMATION OF THE form as V -> -V under space inversion. Because of 

WEYL HAMILTONIAN UNDER SPACE 
INVERSION 



The ingredients which intervene in the Peierls substi- 
tution (|9|), q — > — iV + eA(r), are vectors, which trans- 
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this transformation property and because the general- 
ized Weyl Hamiltonian ([1]) is linear in the momentum, 
we have 



ff s (r,V) = -F s (-r,-V), 



(Bl) 



in space representation, where Hb is given in Eq. (fTTj) . 
We now consider a solution of the Schrodinger equation, 

# B (r,V)V e (r) = «Mr), 

of energy e. Due to the property (|B1[) under space inver- 
sion, one finds that 

ff B (r,V)Ve(-r) = -H B (-r,-V)Ve(-r) = -e^e(-r). 

(B2) 

This means that ip e (—r) is a solution of the Schrodinger 
equation with energy — e, and to each energy e in the 
upper energy band there exists one at — e in the lower one. 
The resulting Landau level energy spectrum is, therefore, 
symmetric around zero energy. 

Notice, however, that this argument is only valid for 
the generalized Weyl Hamiltonian if one neglects the un- 
derlying lattice model for which the Weyl Hamiltonian 
describes the low-energy properties. As mentioned in 
the text, Dirac points occur in pairs at non-zero wave 
vectors, and one thus obtains a twofold valley degener- 
acy £ = ±1. In this case, there is a relative minus sign 
between the two copies of the Weyl Hamiltonian £,Hb- 
Space inversion, therefore, becomes an exact symmetry 
of the model and involves both valleys. 



sent the ladder operators a and as 

ae lv = d x + x, aV^ = -d x + x , 

which act on states described as functions of the real 
variable x. The latter is not necessarily identical to the 
x- or x'-variable in the plane, but should rather be viewed 
as an auxiliary variable that allows for a formal solution 
of the zero-energy problem. 

The Hamiltonian (fTTj) . thus, reads 



H R = 



^J2w x w y f w x (x + d x )e 



h 



{x-d x y? 



WqX 



, (CI) 



and we may search the solutions of the zero-energy eigen- 
value equation, HbtPo(x) — 0, with the ansatz 



This yields the system of linear equations 

w a + e"^(l = 

e^(l+7)" + w (3 = , 
which has a non-zero solution if we choose 



(C2) 



and, thus, 



APPENDIX C: ZERO-ENERGY LANDAU LEVEL 
FOR THE CASE WITH TILT 

In this appendix, we show that there exists a zero en- 
ergy mode for the Hamiltonian Hb (fTTj) . We may repre- 



w e 



1 + 7 



-a . 



Notice that we only obtain bound states, i.e. states that 
decay at large values of x, for wo < 1 (7 > and real), 
which corresponds to the maximal-tilt condition ([5]). 
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